Probing active forces via a fluctuation-dissipation relation 
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We derive a new fluctuation-dissipation relation for non-equilibrium systems with long term memory. We 
show how this relation allows one to access new experimental information regarding active forces in living cells 
that cannot otherwise be accessed. For a silica bead attached to the wall of a living cell, we identify a cross-over 
time between thermally controlled fluctuations and those produced by the active forces. We show that the probe 
position is eventually slaved to the underlying random drive produced by the so-called active forces. 
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Living cells are paradigmatic out of equilibrium systems, 
subjected to the ATP-driven activity of a collection of molec- 
ular motors, whose individual motion cannot easily be disen- 
tangled from thermal fluctuations. These are relatively small 
systems for which fluctuation phenomena are prominent, as 
investigated in recent works [1<h4], with specific focus on 
the active forces spectrum fl] 15) which can be measured via 
micro-rheological devices. Our goal in the present work is to 
show how the cell tissue's random pull-and-push can be in- 
vestigated by means of very recent theoretical advances in the 
field of non-equilibrium statistical mechanics. And indeed, 
from the theory standpoint, in the recent past, much effort 
has been invested in deriving simple generalizations or ex- 
tensions of the celebrated fluctuation-dissipation theorem for 
systems that can be arbitrarily far out of equilibrium ll3ll6 HT3l . 
These efforts have given birth to a flurry of formulas relating 
the response of a system to a small external perturbation to 
some correlation functions, even when the system is not in 
an equilibrium state. Yet, so far none of these formulas has 
been of any predictive power in an actual experimental sys- 
tem. Existing experiments revolving around these theoretical 
advances have been confined to refined and nontrivial confir- 
mations that in some small scale systems such as an optically 
trapped Brownian particle Ifl4ifl6l the various ingredients en- 
tering these extended fluctuation-dissipation relations (EFDR) 
can indeed be measured. Our purpose in this work is to ex- 
ploit an EFDR to access relevant pieces of information on ac- 
tive forces in living cells that cannot otherwise be accessed. 
Given that the dynamics of living cells exhibit strong mem- 
ory effects, we will first have to derive our own version of an 
EFDR adapted to a system with stochastic yet non-Markovian 
dynamics. The latter EFDR and its consequence for the un- 
derstanding of active forces are the central topics of this letter. 
For instance, how relevant is the long term memory in relat- 
ing response and fluctuations out of equilibrium (the answer is 
fully negative in equilibrium)? What new pieces of informa- 
tion on non-equilibrium forces can be learnt from such a rela- 
tionship? Can thermal fluctuations be disentangled from those 
arising from the active forces in a quantitative way? These are 
the questions we wish to address in the present work. 

Our experimental setup consists of a pre-muscular (C2C12) 
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FIG. 1. Schematic view of the bead probe as attached to the cy- 
toskeleton. The active force are the by-product of the ATP-driven 
molecular motors activity. 



cell to which we slave a micron sized silica probe specifically 
bound to the membrane. Its dynamics strongly feels the visco- 
elastic non-Newtonian underlying medium [1|, which gives 
rise to memory (non-Markovian) effects. A schematic view 
of our experimental setup is shown in FIG. [T] We assume 
that the bead dynamics is well described by a visco-elastic 
description. Our bead attached to the cell cortex thus evolves 
under the combined effect of thermal fluctuations £(£) and ac- 
tive forces F a (t). Discarding inertial effects, the position x of 
the bead is governed by the Langevin equation 



m / dt'-f(t-t') 



dx(f) 
dt' 



F a (t)+m 



(1) 



where the Gaussian colored noise £ has correlations cr(r) = 
(£(£)£(£ + t)). We assume that the causal function 7 is re- 
lated to a through a local detailed balance condition |[T7l : 
<t(t) = m/^7(|r|), where /3 = l/k B T is the inverse tempera- 
ture. The active forces are modeled by an isotropic stationary 
random force F a (t), While the presence of the memory ker- 
nel 7(7- ) looks like an innocuous extension of the well-known 
overdamped Langevin equation (for which 7(7-) is to be re- 
placed by a Dirac delta in time), we shall soon see that a num- 
ber of complications emerge. 

We begin with deriving our central theoretical result, 
namely a relationship between three quantities, two of which 



are accessible to measurement. Let R(t) be the response of 
the bead position at time t when applying an external space 
independent force at time 0. Let also C(t) — (x(t)x(0)} be 
the auto-correlation function of the bead's position. Our goal 
is to connect R and C to the statistics of the active force F a . 
Then, we shall show how to exploit our mathematical identity 
to extract useful pieces of information on active forces. Start- 
ing from ([TJ, the probability weight for a given realization of 
the noise is written as P[£\ — exp(— S[£\), where the action 
functional S reads: 

S[$ = \fdti [dt 2 Z(h)a- l (t 1 -t2Mt2) . (2) 

The function er _1 is the deconvolution inverse of the noise 
covariance: f dt<r(t)<r~ l (T — t) = S(t). The next step is 
to rewrite dzl by replacing the noise £(£) with its explicit ex- 
pression which solves (QJ. Then, by considering the effect of 
a small force f(t) applied to the system, one has a variation 
from S to S + SS, where 5S is, up to linear order in /: 



SS = dti \ dt 2 \ dt' 1 [-pa(t 1 -t' 1 )v(t' l )(r- 1 (t 1 -t2) 



x f(t 2 )] + [Fairer- 1 (t! - t 2 )f{h)\ + O(f) . (3) 

By defining the response R to the perturbation / as S(x) (t) = 
J dt'R(t - t')f(t') one finds: 
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P Ids I dra(s - r)cr _1 (s - t)-r(x(r)x(t)) 



dsa 



dr 



(s-?)(F a (s)x(t)) . (4) 



Formula Q neither simplifies in the stationary (non- 
equilibrium) limit nor in the absence of active forces F a be- 
cause of the memory-effects. In the Markov limit, the memory 
kernel is a Dirac delta j(t) = jS(t), which leads to 

R(t-t') = ^{x(t')x(t))-^{F a (t')x(t)) , (5) 
2 at' mj 

a formula that as such first appeared in lfT2l and for which 
the second term in the right hand side has an interpretation 
in terms of dynamical activity. The formula being valid for 
arbitrary forces F a , it would also be valid if F a were an equi- 
librium force field, and in that case the terms in the right 
hand side of |5) add up and yield the standard form of the 
fluctuation-dissipation theorem. Finally, in the absence of an 
external force field F a , no stationary regime can be reached, 
but the formula R{t—t') — | ^ (x(t')x(t)} suggests an effec- 
tive temperature twice that of the bath governs the dynamics, 
a simple interpretation already noted in 1 18]. 

The gist of equation Q is that both R(t — t') and the com- 
panion correlation function (x(t')x(t)) can independently be 



measured, and thus, provided the functions cr _1 , a, 7 en- 
coding memory effects are well-known, access is granted to 
the projection of the active force onto the beads position, 
{F a (t')x(t)). Our only requirement is thus to be able to de- 
termine the functions 7 and a~ l . This can explicitly be car- 
ried out by computing the effect of the external perturbation 
5{x) (t) from (fill. For this it is useful first to define the inverse 
memory kernel 7 -1 such that /7(^)7 _1 (i — r)dt = 5(t). 
With this we can explicitly solve (fTl for x(t) obtaining: 



R(t) 



drj (t) 



(6) 



Before precisely showing how to extract useful informa- 
tion from experimental data, we shall first consider the time 
dependence of our correlation functions. Since our system 
is not in a steady state, all correlation functions a priori de- 
pend on two times. However, since the velocity dx(t)/dt is in 
a steady-state, the mean-square-displacement Ax 2 (t,t') — 
([x(t) — x(t')} 2 ) only depends on t — t'. We can then rewrite: 



dt 



7 {x(t')x(t)) 



1 d Tl ^{t')) + \^x\t-t'). (7) 



2d? 



2dt 



Furthermore, the response function only depends on the time 
difference t — if, so that the correlation function (F a (t')x(t)) 
can also be decomposed as: 



(F a (?)x(t)) = (F a (?)x(t')) + C xF (t - 



(8) 



We can now deduce the active-force-position correlations 
from Q, which we can rewrite as a first relation for the initial 
time t' and a second one for the time difference r = t — t': 



(Fa(t')x(t')) = £ 



P 



dra(t'-r)-(x 2 (r)) , (9) 
dr 



C, ,,( '■) = !/ dta(t~T)j t Ax 2 (t)- I dtR(t)a(r-t) 



(10) 



These formulae allow us to fully determine the time behavior 
of the correlation function (F a (t')x(t)}. Of course the phys- 
ically relevant component of this correlation function is the 
time-translational invariant one, denoted above as C x p . In 
fact the first term in the rhs of (TTOl only expresses the de- 
pendence on the choice of the initial time t' . Here we shall 
restrict ourselves to the case where ? = (the general case 
t' y^ will be discussed elsewhere lfl9l ). Furthermore, the 
second term in the rhs of ( fT0| ) corresponds to the correlation 
function between the position and the thermal noise. 

We shall now apply our result to probe active forces in liv- 
ing cells. The experimental set-up is the same as exposed 
in |@) and consists of a small bead attached to the cortex of 
a muscular cell (see FIG[T]l at room temperature (~ 298K). 

First, we track the free diffusion of the bead on the cell 
cortex. This allow us to measure the mean square displace- 
ment of the probe as a function of time (Ax 2 (i)), shown in 
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FIG. 2. The probe mean square displacement as a function of time 
when there is no external applied force. A definite crossover time is 
found at r ~ 1 s. 



FIG. 3. Average particle position (black) and applied force (red) in 
dimensionless units (d is the particle diameter), as a function of time. 
Data is compared with the best fit used in deconvolving operations to 
restrain noise amplification. 



FIG. [2] It exhibits a two-step growth characterized by an early 
sub-diffusion, up to a 1 s. timescale yet to be interpreted, 
followed by super-diffusion. Second, the response function 
is measured by applying optical tweezers to the bead. The 
tweezer acts on the bead as a harmonic potential of spring 
constant k opt = 120pN//im. After switching our tweezers on 
we quickly translate them by xq — 0.6/im at time t — 0, cre- 
ating a displacement of the bead attached to the cell. In the 
experiment the external force we apply is hence well approx- 
imated by fex P [x(t),t] = -k opt (x(t) ~ X )9(t), assuming 
that the displacement at t = is instantaneous. The distance 
between the bead and the center of the tweezers allows us to 
measure the force applied on the bead. By repeating this ex- 
periment we can measure the average displacement and the 
average force as shown in FIG. [3] The force suddenly in- 
creases at a time of about 0.004s and then slowly decreases 
towards 0, while the bead displacement increases following 
a power law behavior (except in the initial steps of the relax- 
ation), as shown in the Log-Log scale inset. This is compatible 
with the visco-elastic rheology responsible for the subdiffu- 
sivity of the mean square displacement [20-22]. The analysis 
also shows that a sum of two power laws fits extremely well 
the data over all the experimentally accessible range. Equa- 
tion (fTl implies that (a;) and (f exp [x}) are related according 
to S (x(t)) = / dt'R(t - if) (f exp [x(t'),t'}), where R is the 
response to a space independent force, as obtained in (BY 

All our experimental data are discretized in time. We 
first measure trajectories of the bead position when no ex- 
ternal force is applied. We have collected Ng — 39 inde- 
pendent sample trajectories of 250s each with a time step 
of 0.004s (each trajectory consists of N p = 62500 points). 
This yields a set of vectors x\ , where the subscript i stands 
for the time position and superscript j refers to the trajec- 
tory number. We compute the mean square displacement as 
(Ax 2 ) t = Ej,fc(4 " x k J+l ) 2 /N a s {N° - i), from which we 
take the discrete derivative. A second set of measurements is 
taken when we apply the external force on the bead. Then 



we measure both the position x(t) and the external force f(t), 
as shown in FIG. [3] All trajectories are measured for a time 
window of 25s with a time step of 0.004s and ensemble av- 
erages are taken over 39 realizations. We obtain the discrete 
response function Ri by deconvolving the discrete external 
force (/) with 8(x). In order to avoid noise amplification dur- 
ing the deconvolution procedure, we have used the fitted data 
instead of real data. From Ri we can access the triangular ma- 
trix 7« , which is the discrete equivalent of 7 _1 appearing in 
(roll. Then 7^ is obtained by numerically inverting "fZ , and 
<r%j =m/f3('j ij +Jji). 

It is now a matter of simple data manipulation to obtain the 
correlation functions (F a (0)x(t)) and (£(0)x(i)), which are 
plotted in FIG. HI As expected, the thermal noise-position 
correlation function is of the order of k B T and decreases to 
after a time of the order of the second. Besides, the active 
force-position correlation function has a completely different 
behavior. First, we observe that the short time behavior ex- 
hibits small yet significant negative correlation. This would 
mean that on average the particle moves opposite to the active 
force for about the first ms. At later times, as expected the cor- 
relation function turns positive and grows linearly with time. 
This confirms the more intuitive idea that the particle on av- 
erage moves in the same direction of the active forces. How- 
ever we can now quantify this correlation, which turns out to 
be of 10 — 100fc B T. The figure also shows that for short time 
(< Is) the effect of thermal noise prevails, while at large times 
the active force dominates. Again the crossover timescale be- 
tween these two regimes is of the order of 1 s. This seems 
to corroborate a scenario for which the short time behavior 
has equilibrium-like properties, from which the sub-diffusive 
nature of visco-elasticity emerges, while the long time be- 
havior is strongly governed by active non-equilibrium forces. 
Finally we briefly comment on the initial anti-correlation ex- 
hibited by the active force-position correlation function. This 
feature appears for the very first milliseconds of the measure- 
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FIG. 4. Position-force correlation as a function of time difference. 



merits. We also add that in this time range the data fit (as 
shown in FIG. [3]) is less accurate than at later times. This 
anti-correlation could however be explained with the help of 
the depicted scenario. Indeed, if the short time dynamics can 
be seen as a diffusion in a localized potential determined by 
the local structure of the cell cortex, the effect of the active 
forces at these timescales could be a deformation of this po- 
tential, which would result in a systematic tilting in the direc- 
tion opposite to the force. Another useful information that can 
be extracted from the correlation function concerns the power 
spent by the active forces on the bead. This can be estimated 
from the time derivative d/dt (F a (0)x(t)) when t — > 0. Our 
experimental data shows that this power is about 30k B T/s. 
Assuming a ballistic-like behavior, we estimate from FIG. [2] 
a typical velocity of v ~ 5 x 10 -3 /xm/s (compatible with 
0.08/im/s observed in l23l ). This leads to an applied force of 
almost 25 pN, which would mean that a small number (five or 
less) of motors are contributing to the bead motion. 

The equilibrium fluctuation-dissipation theorem relates the 
response of the system to a small perturbation to sponta- 
neous fluctuations in equilibrium. Out-of-equilibrium exten- 
sions involve, besides dissipative aspects, also kinetic aspects 
where the dynamical activity appears within an extra corre- 
lation function with the physical observable under scrutiny. 
For overdamped Langevin dynamics, and even for dynamics 
with strong memory effects, the latter is directly related to the 
non-equilibrium forces driving the system out of equilibrium. 
Understanding the properties of these forces in the dynamics 
of living cells is a problem of its own. With this new theo- 
retical tool we have shown how to access some of the prop- 
erties of these forces. We have been able to quantify the time 
scale at which active forces driven fluctuations win over ther- 
mal ones. Another quantitative spin-off is an estimate for the 
power dissipated by the active forces into the system, which 
turns out three orders of magnitude smaller than the chemi- 
cal power injected into the underlying motors by the ATPase 
steps. Another conclusion that we have arrived at is the strong 



slaving of the probe motion to the active forces, which are co- 
herently followed. We have thus demonstrated that extracting 
useful information is indeed possible, and this paves the way 
to a number of improvements or generalizations. An immedi- 
ate modification of the experimental setup would involve the 
use of an optical trap to render the unperturbed state station- 
ary, and then pulling on the particle by displacing the trap. 
Making the system stationary would eliminate time drifts and 
would simplify the subsequent theoretical analysis by allow- 
ing to focus solely on active forces. Other simple improve- 
ments of the present work would, for instance, entail tracking 
higher moments of the bead positions. This could be done 
by attempting to vary, in the latter stationary setup, the trap's 
stiffness, instead of the position of its minimum. Even if the 
rheological meaning of the response function would then be 
lost, other relevant information concerning fluctuations would 
then be gained. We leave these projects for future work. 

We acknowledge several useful discussions with Julien 
Tailleur and Damien Robert. 
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